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Using a novel finite size scaling Monte Carlo method, we 
calculate the four, six and eight point renormalized coupling 
constants defined at zero momentum in the symmetric phase 
of the three dimensional Ising system. The results of the 2D 
Ising system that were directly measured are also reported. 
Our values of the six and eight point coupling constants are 
significantly different from those obtained from other meth- 
ods. 



I. INTRODUCTION 



The Hamiltonian of the Ising ferromagnet is given by 



H 



(1) 



<hJ> 



where the Ising spin at site i, Si, can take either 1 or 
-1 and the sum is over all the nearest neighbors of the 
lattice. It is well known that the critical behavior of the 
D dimensional Ising model can be described by the D di- 
mensional Euclidean scalar field theory the Hamiltonian 
of which is given by 



H = 



d D 3 



, (2) 



where too and go are respectively the bare mass and cou- 
pling constant defined in the absence of critical fluctu- 
ations of the fields. Near the criticality mo is a linear 
measure of temperature, and we denote by mo c the value 
for which the theory becomes critical. As the fluctua- 
tions become strong renormalizations of the mass, cou- 
pling constant and fields are necessary, and the long dis- 
tance behavior is no longer described by the bare poten- 
tial but by the effective potential which is generally of 
more complicated functional form than the bare one. In 
statistical physics the effective potential represents the 
free energy density as a function of order parameter (ex- 
pectation value of the renormalized field) , and is used to 
determine the equation of state. 

After small-renormalized-field expansion of the effec- 
tive potential, its coefficients are directly related to the 
renormalized coupling constants (RCCs) defined at zero- 
momentum. In terms of the expectation value of the 
renormalized field, (fR, the effective potential in the 3D 
symmetric phase may be written as 

Veff(<pn) = ~mfj¥>| + im fl ^Vfl + 



1 1 (8) 

6\ 9r tPR+ 8\m R (PR + 



(3) 



where wir and g^ represent respectively the renormal- 
ized mass (inverse correlation length £) and the (dimen- 
sionless) A-point RCC defined at zero momentum. 

The formal expression of g^ can be obtained by cal- 
culating the A— th derivative of the effective potential 
with respect to the average value of the renormalized 
field. Necessary elements for the calculations are the 
well-known relations 



dV f 



eff 



dp 



J, ip 



1 dW[J] 
V dJ ' 



and 



(4) 



Here ip, V, W, and Z v are respectively the expectation 
value of the bare field 4>, the volume of the system, the 
generating functional for connected Green functions in 
the presence of the external field J, and the field strength 
renormalization factor given by Z v = x m R with \ denot- 
ing magnetic susceptibility. The expressions in 3D are 



(4) 

9r 

(6) 

9r 

(8) 

9r 



-(Z^/m R ) W 2 - 4 W A 



-Z% Wi 6 [W 6 - 
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lowf W 2 X \ 
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(5) 
(6) 

(7) 



where Wn is the Fourier transformed A-point connected 
Green function at zero momentum. In 2D ip R is dimen- 
sionless, so that the expressions of the g R N ^ are given by 
the above equations with additional division of the m R 
term. 

Each RCC has a universal zon-zero finite value as the 
theory becomes critical, provided the hyperscaling rela- 
tion holds. To be more specific, we take thermodynamic 
limit first and then take the limit of the thermodynamic 
correlation length becoming divergent (that is, the limit 
too — > TOo c ). The corresponding value in this limit may 
be called as critical RCC. The determination of the accu- 
rate value of the critical four-point RCC (g R ) is partic- 
ularly important since all other universal quantities can 
be given in terms of it. 

In this work we are mainly concerned with the Monte 
Carlo calculation of the critical RCCs in the symmetric 
phase of the three dimensional (3D) Ising model. (From 
here on the notation of RCC will be used without the sub- 
script R when it is used in context of size dependence.) 
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This has been a subject of many studies ^including a vari- 
ety of quantum field thettrgtic approachscTu, high temper- 
ature series exuapsionsolD, exact renormalization group 
flow techniquesQ~tj, and Monte Carlo simulationsOEJ'EI 
It is known that standard direct Monte Carlo measure- 
ments of the thermodynapuc critical RCCs suffer from 
enormous statistical noiseO. Other Monte Carlo study 
is based on the finding of the probability distribution 
of the order .parameter in the external field at various 
temperaturealj. . We employ a novel finite size scaling 
(FSS) techniqueEaEJ combined with single cluster flip- 
ping Monte Carlo algorithmic. Our preliminary results 
were reported in Ref.( |l5|). Here we report the result of 
our updated Monte Carlo measurements that has been 
increased more than one magnitude of order. We also 
report our results measured at criticality that is useful to 
determine the correction to scaling in the RCCs as well 
as results in the 2D system. 

II. THE MONTE CARLO METHOD 

For the calculations of the RCC one needs to calcu- 
late £, X; an d W2N ■ With periodic boundary condition 
imposed on the lattice, x an d W2N can be expressed in 
terms of the expectation values of various powers of the 
sum of the spin over all the lattice site S = Y^j F° r 
example, 

X =< S 2 > I L D , W 4 = (< S 4 > -3 < S 2 > 2 ) /L D , 

(8) 

where L is the linear size of the lattice. £ can be very ac- 
curately-^determined using the stapdard second moment 
formulate. The cluster algorithmsE-3 of the Monte Carlo 
simulation have been extremely efficient for many prob- 
lems of critical phenomena. It is now an easy task to 
obtain Monte Carlo data of typical physical quantities 
such as x an d £ with relative statistical errors less than 
0.1 percent at a temperature arbitrarily close to critical- 
ity. Nevertheless, it turns out that accurate Monte Carlo 
measurement of higher order RCC is problematic in some 
cases. The problem stems basically from the fact that the 
higher order RCC is given as a multiplication of a huge 
number with a tiny one. The former is given typically 

by some power of for instance ( ^-J for the g^' of 

the 3D system, and the latter comes from the combina- 
tion of the W2N that turns out to be extremely sensitive 
to statistical noises. The noises increase rapidly with in- 
creasing temperature and increasing value of 

To illustrate we measured various physical quanti- 
ties at an arbitrary (inverse) temperature in the scaling 
regime of the 2D Ising model with increasing linear size 
of the lattice L. Table(||) clearly shows that the statisti- 
cal error in the measured values of g^ N ^ increases with 
N for any given value of L. For L=80, for example, the 



relative statistical errors of £l , XL, 9 l ■, 9l j 9l are re ~ 
spectively 0.08, 0.1, 0.7, 3.1, 5.3 percents. We observe 

(8) 

clearly that all the variables except for g K L ' are mono- 
tonically increasing functions of L, up to L=80 where 
L/£l — 6.7. For further increasing L the values do not 
vary within the statistical errors. We recall that the well- 
known exact thermodynamic value of the £ at this (3 is 
11.9055..., which manifests itself at L=80 already within 
the statistical errors. The size independent value thus 
corresponds to the thermodynamic value. It is also very 
clear that the size dependence becomes rapidly weaker 
with increasing L for all the variables we consider here. 
This is expected to be the case for any physical variable 
having well-defined thermodynamic value. For the g L , 
however, due to the large statistical noises is it hard to 
draw a definite conclusion on the size dependence. Never- 
theless, it is very natural to expect that g L would follow 

similar size dependence as g^ and g^ • 

To compute the critical RCC it is required to obtain 
thermodynamic values as (3 — ► (3 C . For the case of the 2D 
Ising model, we will not follow that procedure. Instead, 
we will report the results of measurements at criticality 
and try to find out possible effect of the correction to 
scaling. The illustration in the case of the 2D Ising model 
will be useful in connection with the study of the 3D case. 
Note that the relation in the scaling regime 

ff g JV) (t)~t D "- aA ^ ) (9) 

translates into 

g^ N \t = 0)~L^- 2A ^+^ v , (10) 

at criticality. Hence the hyperscaling relation Di> — 
2A2N + 7 = implies the invariance of g^ (t — 0) with 
respect to L. Conversely, from the (weak) dependence of 
of gf N \t = 0) can one infer the correction to scaling in 
the scaling regime. Computing RCC by taking the limit 
t — > first and then take the limit L —> 00 is equiva- 
lent to the computation of g^ 2N \t = 0). It should be 

stressed that g^ N \t = 0) basically represents the appo- 
site of thermodynamic limit for a non-zero t — ► (that is 
L — > limit for £ becomes huge but not infinite) and-has 
nothing to do with the lower-bound of critical RCCEi. 

We also report so called 6-th and 8-th order cumulant 
ratios at criticality that are denoted by uf^ and uf^ 
respectively. These are part of We and W%, and have al- 
ready been available for the 2D Ising medel based on um- 
brella sampling method of the Monte Carlo simulatior£3. 
Since our code is generic in dimension, the comparison 
with the result from different Monte Carlo method for 
these higher order cumulant ratios having severe statisti- 
cal noises may be useful to check against any unexpected 
possible systematic errors in our Monte Carlo measure- 
ments. For example, is expressed by 

Uf ) = (< S 8 > -28 < S 6 >< S 2 > -35 < S 4 > 2 + 

420 < S 4 >< S 2 > 2 -630 < S 2 > 4 )/ < S 2 > 4 (11) 
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The results are summerized in Tablc(gl|). Our results 
of t/f N) for N=2,3, and 4 are respectively 1.8318(5), 
13.93(1), and -226.0(1), which may be compared with 
1.834, 13.96, and -226.6 reported in Ref. ||. Currently 
accepted value of is 1.832(1). We observe that 

g^ N \t — 0) has no L dependence at least for L > 40, 
confirming the hyperscaling and no significant correction 
to scaling. We therefore lead to the conclusion that the 
thermodynamic values of RCC measured at (3 = 0.420 
may well be regarded as the critical values. 

Regarding size dependence we observe similar feature 
in the 3D Ising model as in the 2D case. Our data are 
summerized in Table Q. For the (3 = 0.217 and 0.220 
we observe no size dependence (within the statistical er- 
rors) for the £ and \ beyond L=32 and 60 respectively. 
This roughly corresponds to L/£l ~ 5.5. We expect that 
this is the case even for the RCCs, although it was almost 
impossible to get precise measurements of the thermody- 
namic values of gf^ and for this value of L/£l. For 
example, for (3 = 0.217 and L = 32 we generated about 
10 9 single cluster sweeps, but error bars are larger than 
the mean value for the g^ . The values of gf^ for L > 32 
also seem to be unreliable in view of the generic feature 
of the weaker size dependence with larger L. 

It thus appears that Monte Carlo computation of the 
critical RCCs relying on direct brute force measurements 
is prohibitively difficult. In order to overcome the diffi- 
culty closer, to T Cl we make use of a new FSS function 
Q A {x(L, t))UH, defined by the expression 

A L (t) = A(t)Q A (x(L,t)), x{L,t)=t L (t)/L. (12) 

Here A^(t) represents the quantity A measured on a finite 
lattice of linear size L at a reduced temperature t, with its 
corresponding thermodynamic value A(t). What Eq.jl^) 
states is that the size dependence of a physical quantity A 
is given as a function of the scaling variable x. As a result, 
the ratio of L to £l beyond which thermodynamic limit 
is reached is independent of the temperature, that turned 
out to be approximately 5.5 for the 3D Ising model. 

The FSS technique is especially useful for our pur- 
pose, because it enables us to extract accurate thermo- 
dynamic values based on the Monte Carlo measurements 
with much smaller lattices. We just outline the single 
step FSS extrapolation technique used in this work. For 
a detailed explanation, we refer the readers to Ref. [l3 14 



1. For a certain £o ? measure Ai,(ta) and x(L,to) = 
£l(^o)/£ for increasing L. 

2. Determine the thermodynamic value at the tem- 
perature A(to) by measuring At (to) which is L in- 
dependent. 

3. Fit (x(L,to), Al(4q) I 'A(to)) data to a suitable func- 
tional form. In this work we used the ansatz, 



4. For any other t, choose a suitable L, measure the 
value of x(L,t) = £l/L and Ai J {t) 1 and interpolate 
Q(x(L,t)). 

5. Extract A(t) by plugging A^it) and Q(x(L,t)) into 
Eq.@. 

The smallest value of L we considered for the FSS method 
is 20. 



III. RESULT AND DISCUSSION 

Our choice of j3o is /3q = 0.220. We infer from the 
generic feature of the size dependence observed for the 
3D Ising model that thermodynamic limit of the RCCs 
are reached for L/£ L > 5.5. We thus get g^ = 24.5(4) 

and g^ = 1983(148) at this (3. Fitting the data to the 
ansatz Eq.(|l3|), we get 

d = 2.338, c 2 = -15.768, c 3 = 19.770, c 4 = -5.123 
ci = 3.225, c 2 = -25.046, c 3 = 42.485, c 4 = -23.177 

for the scaling function Q(x) of g^ and g^ respectively. 
Using the scaling function we calculated the thermody- 
namic values of the 4-point and 6-point RCCs for all the 
L from 36 to 80 in the Table|v|). The result at each L 
is in reasonably good agreement. For example, we get 
g ( ^= 24.3(1), 24.4(4), 24.1(1), 23.9(1), 23.9(1), 24.1(3), 
and 23.9(2) for each L from the L=36 through the L=80 

in the table. For g^ ] we get 1919(11), 1939(20), 1939(21), 
1915(10), 1906(17), 1917(70), and 1897(49). The invari- 
ance of the thermodynamic RCC with respect to the 
choice of L is a numerical proof of the FSS for the vari- 
ables (see Figure(Q) and (||)). We usually extracted the 
thermodynamic value for several different choices of L 
for a given temperature, and took the average. Our net 
results from @ = 0.217 to (3 = 0.2213 are found in Ta- 
blc(|vi|). It is observed that both g^' and <?jj tend to 
decrease mildly as — > C . In this work we assume the 
widely accepted correction to scaling exponent 9 ~ 0.5 
and /3 C = 0.221654^3. By fitting our data in Table® to 



9r W=9r '{l + a 2N t v ") 1 
we obtain the critical RCC which reads 



-(4) 

9r 



23.6(2) 



= 1879(50) 



(14) 

(15) 
(16) 



Q(x) = 1 + c\x + C2X 2 + c 3 x 3 + C4X 4 



(13) 



Our results of g^ N \t = 0) are found in Table©. 

It is observed that both of and £,l/L have ten- 

dency of very mild decreasing with increasing L. All the 
g^ N ^ (t = 0) show remarkable invariance with respect to 
increasing L at least for L > 30. In other words, they 
do not show the effect of correction to scaling observed 
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in the scaling regime. There may be a few possible inter- 
pretations for the discrepancy. First let us remind that 
Eq.(|l4"|) translates into 



0) 



0){l + b 2N L- 



J ), («>0), 

(17) 



and that the coefficient b^N may happen to be very small. 
The second possibility is that currently accepted (3 C might 
be slightly underestimated: At exact criticality it is ex- 
pected that (,l/L remains an exact constant irrespective 
of the presence of the correction, which seems not to be 
the case at /? = 0.221654. 

Our current result at criticality may be summerized as 
follows. 



g (4 \t = 0) = 5.36(2) 



9 

g {6) (t = 0) = 154.6(1.1) 
<? (8) (t = 0) = 1.035(10) x 10 4 . 



(18) 
(19) 
(20) 



It is impossible to apply the FSS method for the g R 
due to its large error bars. Nevertheless, it is highly 



likely that at least g R ' is of order 10 5 . (see Table(IH).) 

Also we have reliable data at /? = 0.2213 showing g L ~ 
1.13(3) x 10 5 already for L/£ L ~ 2.7. Furthermore it is 

observed that the ratio of critical RCCs to g^ N ^ (t = 0) 
tends to increase with N; The ratio is roughly 4.4 and 
12.2 for N=2 and 3 respectively for the 3D. We also note 
that L dependence of g * {t — 0) is quanitatively the 
same as that of the other RCCs, indicating that the cor- 
rection to scaling for g^' (t) is as mild as the other RCCs. 
All the evidence almost undoubtedly point out that the 

(8) 

critical g R is much larger than those estimated by other 
methods. Our crude estimate of g R reads 

3^ 8) a 1.4(3) x 10 5 . (21) 

Our value of g^ R is in reasonable agreement with other 
estimates. The agreement is especially i-good with the 
result from high-temperature expansionO and the field 
theoretic treatment^. However, the agreement becomes 
worse as N increases: The results of g R from previous 
studies are within the range 860 < g R ^ < 1515. The 

closest result to ours, g^ ~ r^l^, obtained from the 
previous Monte Carlo methodEJ is remarkable because 
it is obtained not by the measurement of RCC but by 
the Monte Carlo measurement of the probability distri- 
bution of order parameter assuming the thermodynamic 
limit for L/£z, ~ 4. We conjecture that this result would 
agree much better with our estimate if the same thermo- 
dynamic limit as in this work would have been taken. The 

values of g R from previous studies ranges from 2.9 x 10 4 
to 3.5 x 10 4 , which is at least four times -smaller than 
our estimate. Although it was pointed oura that longer 
series terms are necessary for more accurate estimate of 



higher order critical RCC from the high temperature se- 
ries expansion method, it remains mysterious why previ- 
ous studies based on different methods give rise to rea- 
sonably consistent results among them. Nevertheless, in 
view of limitations of almost all the methods used to 
study this subject, it is fair to say that the issue is still 
open to further studies. 

Our results of the critical RCCs for 2D Ising model 
obtained assuming negligibly small correction to scaling 
read 



The only result on the 
is the e expansion stud; 



14.7(2) (22) 
850(25) (23) 
8.9(5) x 10 4 . (24) 

h order RCCs we are aware of 



~(6) 
9r 



794, 



-(8) 
9r 



.2(2) x 10 4 



which agree reasonably well with ours. 

The results at criticality are summerized as 



~9 (4} (t = 0) 

~9 (6) (t- 
~9 {8) (t = 0) 



2.239(7) 
0) = 29.34(20) 
947(10) 



(25) 



(26) 
(27) 
(28) 



We wish to thank Maxim Tsypin for many communi- 
cations. 



TABLE I. Size dependence of 
at /3 = 0.420 up to L=100 for the 



various physical quantities 
2D Ising system. 



L 


Zl 


XL 


(4) 


(6) 

9l 


(8) 


20 


9.664(8) 


116.1(1) 


6.29(1) 


218.46(1.32) 


17871(140) 


30 


10.977(9) 


162.5(1) 


8.95(1) 


414.8 (2.3) 


43532(313) 


40 


11.53(1) 


186.5(2) 


11.26(1) 


615.1(4.0) 


71443(615) 


50 


11.78(1) 


197.4(2) 


12.83(3) 


740.2(5.5) 


85462(879) 


60 


11.84(1) 


201.3(2) 


13.79(4) 


794.8(6.2) 


87919(1185) 


70 


11.88(1) 


203.1(2) 


14.20(6) 


825.1(10.0) 


89824(2315) 


80 


11.91(1) 


203.9(2) 


14.77(10) 


850.4(26.1) 


89521(4765) 


90 


11.91(2) 


204.3(2) 


14.85(12) 


858.8(24.2) 


90722(7451) 


100 


11.90(2) 


204.4(3) 


14.60(16) 


846(38) 


88165(9575) 



TABLE II. Size dependence of various physical quantities 
at criticality (3 C = In (s/2+ 1) /2 up to L=100 for the 2D Ising 
system. 



U 



(4) 



u 



(f.) 



V 



is, 



(4) 

9l 



(6) 

91 



(8) 

91 



20 1.8324(6) 13.94(1) 

40 1.8321(6) 13.94(1) 

60 1.8317(5) 13.93(1) 

80 1.8318(5) 13.93(1) 

100 1.8316(6) 13.93(1) 



-226.1(2) 2.227(7) 29.03(18) 933(9) 

-226.1(2) 2.237(8) 29.28(21) 945(10) 

-225.9(1) 2.241(7) 29.37(18) 949(9) 

-226.0(1) 2.240(6) 29.36(16) 948(8) 

-226.0(1) 2.239(7) 29.33(18) 947(9) 
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TABLE III. Size dependence of the various physical quan- 
tities at p = 0.217 (the upper part) and f3 = 0.220 (the lower 

(8) 

part). Note that the g L data for > 5.7 become un- 

reliable due to huge error bars and in view of weaker size 
dependence for larger L. This appears to be case even to the 
gf ] for L = 36. 



L 




XL 


( 4 ) 

91 


(6) 

91 


3i 8) x 10- 4 


8 


3.93(0) 


59.49(6) 


9.86(2) 


492.5(1.9) 


5.45(3) 


12 


4.85(1) 


94.16(12) 


13.59(3) 


852(4) 


10.8(1) 


16 


5.299(1) 


114.54(5) 


17.71(2) 


1277(3) 


15.77(9) 


20 


5.488(3) 


124.20(7) 


21.07(5) 


1562(9) 


15.5(3) 


24 


5.573(2) 


128.40(5) 


23.16(8) 


1735(22) 


15.5(9) 


28 


5.605(4) 


130.06(11) 


24.72(11) 


1887(53) 


14.6(3.8) 


32 


5.622(2) 


130.87(5) 


24.88(14) 


2190(86) 


13.3(16.9) 


36 


5.619(19) 


130.96(14) 


25.64(16) 


2860(191) 


83.3(29.1) 


16 


7.85(2) 


228.8(7) 


9.52(5) 


452.1(4.6) 




20 


8.85(2) 


298.1(7) 


11.43(5) 


625.6(5.1) 




24 


9.56(2) 


351.8(1.1) 


13.44(7) 


814.7(8.2) 




30 


10.20(3) 


407.1(1.1) 


16.5(1) 


1123(16) 




36 


10.56(2) 


439.2(1.2) 


19.2(1) 


1372(27) 




40 


10.68(3) 


455.2(1.5) 


21.4(2) 


1501(56) 




50 


10.83(3) 


467.9(1.3) 


23.7(3) 


1795(109) 




60 


10.89(3) 


472.3(1.8) 


24.5(4) 


1983(148) 




70 


10.90(3) 


473.0(1.1) 


25.8(1.6) 






TABLE IV. 


Size depenc 


lence of the 


: various physical quan- 


tities at P = 0.221 up to L/£ L ~ 4.11 






L 


iL 


XL 


(4) 

9l 


(6) 

9l 


»i 8) x 10- 4 


20 


10.95(2) 


426.7(1.4) 


7.6(1) 


301(3) 


2.67(4) 


28 


13.93(6) 


703.4(4.0) 


9.2(1) 


423(8) 


4.3(1) 


36 


15.98(6) 


944.2(4.7) 


11.2(1) 


600(10) 


6.6(2) 


40 


16.74(5) 


1045.5(4.3) 


12.4(1) 


710(13) 


8.2(2) 


48 


17.91(3) 


1210(3) 


14.59(6) 


920(8) 


10.8(2) 


56 


18.60(1) 


1316(1) 


16.74(5) 


1116(7) 


12.2(2) 


64 


19.03(2) 


1386(2) 


18.66(7) 


1297(13) 


13.4(4) 


72 


19.28(5) 


1426(6) 


20.47(29) 


1465(57) 


13.3(2.6) 


80 


19.46(2) 


1458(2) 


21.36(16) 


1580(82) 





TABLE V. Size dependence of various physical quantities 
at criticality /3 C = 0.221654 up to L=80 for the 3D Ising sys- 
tem. Here represents fourth-order cumulant ratio. 

L W Z XL gj 4) W fffxlO- 4 

20 1.418(1) 12.80(1) 545.1(5) 5.40(1) 157.4(8) 1.07(1) 

30 1.409(1) 19.20(2) 1215(1) 5.37(1) 155.3(7) 1.04(1) 

40 1.408(2) 25.65(4) 2148(4) 5.34(2) 153.5(1.1) 1.03(1) 

50 1.403(1) 31.98(4) 3313(4) 5.36(2) 154.4(1.0) 1.03(1) 

60 1.401(2) 38.33(6) 4742(10) 5.37(2) 155.1(1.1) 1.04(1) 

70 1.399(3) 44.69(8) 6418(19) 5.38(2) 155.2(1.2) 1.04(1) 

80 1.398(2) 51.13(9) 8351(19) 5.36(2) 154.1(1.3) 1.03(1) 
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TABLE VI. Thermodynamic values of the four and six point RCC extracted by the FSS technique for some temperatures over 
0.217 < /3 < 0.2213. Here the quoted errors are obtained ignoring the statistical errors in the estimate of the thermodynamic 
values at j3 — 0.220. Crudely speaking, increases in the estimate of beta = 0.220 would lead to more or less similar amount of 
increase in the estimate for other temperatures. 
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FIG. 1. Q g 4 calculated for ,3=0.220 and 0.221. The two 
data set collapse unto a single universal curve, showing a nu- 
merical evidence for the FSS Eqjl^ for . 



FIG. 2. The same as in Figure(|) but for Q g 
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